Quantum Chaos, Delocalization, and Entanglement in Disordered Heisenberg Models 
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We investigate disordered one- and two-dimensional Heisenberg spin lattices across a transition 
from integrability to quantum chaos from both a statistical many-body and a quantum-information 
perspective. Special emphasis is devoted to quantitatively exploring the interplay between eigen- 
vector statistics, delocalization, and entanglement in the presence of nontrivial symmetries. The 
implications of basis dependence of state delocalization indicators (such as the number of principal 
components) is addressed, and a measure of relative delocalization is proposed in order to robustly 
characterize the onset of chaos in the presence of disorder. Both standard multipartite and gen- 
eralized entanglement are investigated in a wide parameter regime by using a family of spin- and 
fermion- purity measures, their dependence on delocalization and on energy spectrum statistics be- 
ing examined. A distinctive correlation between entanglement, delocalizatton, and mtegrabiUty is 
uncovered, which may be generic to systems described by the two-body random ensemble and may 
point to a new diagnostic tool for quantum chaos. Analytical estimates for typical entanglement of 
random pure states restricted to a proper subspace of the full Hilbert space are also established and 
compared with random matrix theory predictions. 

PACS numbers: 03.67.Mn, 03.67.Lx, 05.45.Mt, 24.10.Cn 



I. INTRODUCTION 

The emergence of non-integrable behavior in quantum 
mechanics is a fascinating and widespread phenomenon 
which is largely responsible for the "complexity" intrinsic 
to the physical and mathematical description of interact- 
ing many-body quantum systems. A most striking con- 
sequence is the existence of distinctive quantum- chaotic 
properties for dynamical systems which may lack a clear 
classical limit. The characterization of such quantum 
chaos signatures has a long history, pioneered by Wigner 
in his effort to quantitatively model complex nuclei 
and eventually culminating in statistical approaches to 
complex quantum systems based on so-called Random 
Matrix Theory (RMT) 

Recent years have witnessed a renewed interest to- 
ward qualitatively reassessing and quantitatively explor- 
ing many-body quantum complexity and quantum chaos 
implications in the light of Quantum Information Sci- 
ence (QIS) [EH]- On one hand, a deeper understanding 
of quantum chaos and its implications is a prerequisite 
for identifying potentially harmful consequences as well 
as beneficial uses of chaos in information-processing de- 
vices: While the possibility that disorder may destabi- 
lize quantum computation through a "chaotic melting" 
0, Q calls for careful hardware design and error control, 
chaotic evolutions tend to naturally generate effectively 
random states which are a resource for a variety of QIS 
protocols [9]. On the other hand, QIS provides additional 
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tools for describing complexity of states and evolutions, 
which are proving useful in a variety of settings at the 
interface with condensed-matter and statistical physics. 
The notion of entanglement, in particular - as captur- 
ing distinctively quantum correlations which admit no 
local classical interpretation - plays a central role to 
this end. Among the most notable developments to date, 
entanglement theory has allowed to considerably deepen 
the understanding of quantum critical phenomena - ap- 
propriate entanglement measures serving as an "order 
parameter" for detecting and classifying quantum phase 
transitions in matter ^11] - and to devise enhanced com- 
putational methods for both static andjime-dependent 
properties of quantum lattice systems [l2l |. 

In this context, uncovering the relationship between 
various aspects of complexity in quantum states and evo- 
lutions, quantum chaos, and entanglement is both a nat- 
ural and fundamental challenge, which is spurring sig- 
nificant activity in the field, see e.g. [H, TT3 . [Tsl . [Ta . 
[13, [H, [H, [13, [m, [13, [H for representative contribu- 
tions. In particular, specific questions to be answered in- 
clude the following: What entanglement properties best 
capture the structural change occurring in typical many- 
body eigenstates across a transition to quantum chaos, 
and how well do such properties reflect the complexity 
of chaotic eigenstates? To what extent does the amount 
of entanglement relate to the amount of underlying state 
delocalization? Perhaps most importantly, can entangle- 
ment theory suggest new signatures of quantum chaos? 

Our goal in this work is to take a step toward answer- 
ing some of the above questions, by seeking an in-depth 
characterization of entanglement properties of the sta- 
tionary states (eigenvectors) of non-integrable as opposed 
to integrable many-body Hamiltonians, in relation to the 
behavior of traditional complexity indicators related to 
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RMT spectral statistics and delocalization measures. In 
this respect, our analysis shares some motivation with 
earlier studies of entanglement across a transition to non- 
integrability in a class of one-dimensional Harper Hamil- 
tonians by Lakshminarayan et al. f24^ . and spin- 1/2 lat- 
tice systems by Santos et al. [IH], and Mejia-Monasterio 
et al. p6j where, however, primary emphasis to pair- 
wise and bipartite entanglement is given. Only recently 
has genuine multipartite entanglement started to be ad- 
dressed, notably in the one-dimensional Ising model with 
a tilted magnetic field [lll,[23|. 

Here, we focus on a representative class of disordered 
Heisenberg models, which have received limited attention 
to date in spite of their prominent role in condensed- 
matter physics as well as in exchange-based circuit-model 
[H, m, [lo, [U, m, m, Hl] and cluster-state ^ quan- 
tum computing architectures. As a further distinctive 
feature of our work, the notion of Generalized Entangle- 
ment (GE), introduced by Barnum et al. in [H, [s^l, is 
exploited to both obtain a unified approach to standard 
(qubit-based) multipartite entanglement - quantified by 
a family of coarse-grained (spin) purity measures - and 
to construct GE measures directly probing correlations 
in different (fermionic) operator languages. From this 
point of view, the present study further validates the 
usefulness of GE for broadly characterizing complexity 
in quantum systems, as recently demonstrated in appli- 
cations to ground-state quantum phase transitions [38j . 
chaotic quantum maps [39|, and efficient solvability of 
Lie-algebraic models pO| . 

The content of the paper is organized as follows. We 
begin in Sec. II by recalling the essential RMT back- 
ground, along with well-established spectral signatures 
of quantum chaos and measures of pure-state delocal- 
ization. Sec. HI introduces the relevant class of one- 
dimensional (ID) and two-dimensional (2D) Heisenberg 
models, laying out the static disorder settings under ex- 
amination (associated with randomness in the one-body 
energies, the two-body interactions, or both), and dis- 
cussing the symmetries associated with different param- 
eter regimes. A thorough characterization of integrabil- 
ity properties is obtained in Sec. IV, with the twofold 
objective of distinguishing between delocalized-localized 
regions and chaotic- integrable ones. The implications of 
the dependence of delocalization upon the basis choice 
are elucidated, and a measure of relative delocalization 
is introduced to gain additional insight into the proper- 
ties of eigenvectors in disorder-induced chaotic regimes. 
After a brief account on entanglement and GE measures 
in Sec. V, we present in Sec. VI a detailed analysis of 
entanglement as a function of disorder strength, energy, 
as well as delocalization properties. As a main emerg- 
ing feature, a strong and persistent correlation between 
multipartite entanglement and delocalization is found in 
non-integrable regimes, which is consistent with indepen- 
dent evidence in [27| and may suggest a novel signature of 
quantum chaos. While most results are numerical, ana- 
lytical estimates for multipartite entanglement of random 



pure states localized to a subspace are also obtained and 
compared to RMT predictions. Throughout the paper, 
special care is devoted to contrast properties which are 
general to disordered spin- 1/2 lattices to those which are 
specific to the Heisenberg interaction. The paper con- 
cludes with a summary and outlook in Sec. VII, followed 
by an Appendix which collects technical derivations. 

II. SIGNATURES OF QUANTUM CHAOS 

It is well established that quantized versions of classi- 
cally integrable and fully chaotic systems can be distin- 
guished by their quantum energy level statistics P, |3, 0, 
EH, 113, E^l . Of particular interest is the distribution of 
energy level spacings, P{s), where s is the spacing be- 
tween neighboring energy levels after the spectrum has 
been appropriately unfolded, so that the density of states 
is everywhere equal to one. Integrable systems typically 
exhibit a Poisson distribution, 

Ppis) = exp(~s), s e R+, (1) 

whereas chaotic systems have an energy level spacing 
distribution predicted by RMT. Within RMT, an ex- 
act description of a complex physical system is replaced 
by a statistical description based on ensembles of ran- 
dom matrices which share the same fundamental sym- 
metry properties as the original system Hamiltonian. In 
particular, for the wide class of systems exhibiting time 
reversal invariance, the appropriate ensemble is the so- 
called Gaussian Orthogonal Ensemble (GOE), whose en- 
ergy level spacing distribution is closely approximated by 
the Wigner-Dyson distribution, 

TTS ( H S^ \ 

Pwd(s) = Y exp ( ) , s e R+. (2) 

The eigenvectors of fully chaotic systems may also be 
described statistically using RMT. If the system of inter- 
est obeys time-reversal symmetry, the eigenvector com- 
ponents tend to follow a Gaussian distribution, resulting 
from the invariance of the GOE under arbitrary orthog- 
onal transformations pj| . 

It is important to appreciate that systems such as in- 
teracting lattices of spin-1/2 particles have no obvious 
classical limit, thus the question of whether or not they 
may exhibit chaos must be posed and answered from an 
entirely quantum-mechanical perspective. Because the 
standard defining features of classical chaos (phase-space 
ergodicity and exponential divergence of neighboring tra- 
jectories) have no direct meaning for Hamiltonian quan- 
tum dynamical systems, neither does an unambiguously 
established framework exist for consistently defining in- 
tegrability in quantum settings [H, , such an issue is 
not straightforward and still largely open. In the present 
context, we shall use the term "quantum chaos" in an 
operational sense, to simply mean the presence of RMT 
energy level statistics. 
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In order to examine the transition from integrability 
to chaos, the foUowing picture is employed [26l. l4ll. [46| . 
Let Hq be an integrable Hamiltonian, that is, one for 
which the eigenvalues and eigenvectors may be deter- 
mined analytically [i^ H^]. Now consider the effect of 
an integrability-breaking perturbation, H' , so that the 
total Hamiltonian becomes 



(3) 



For sufficiently small values of the parameter A, the 
eigenvalues and eigenvectors of H are adequately de- 
scribed by perturbation theory. As A increases to 
the point where the interaction strength between states 
which are directly coupled by H' is equal to their un- 
perturbed energy difference 0, Q, perturbation the- 
ory breaks down, and a crossover from a Poisson to a 
Wigner-Dyson distribution occurs. In parallel with such 
a crossover in energy level statistics, the eigenvectors of 
H become increasingly delocalized across the eigenstates 
of Hq. As A increases further, delocalization typically 
continues until the component distribution of the eigen- 
vectors becomes that of GOE random states. The exact 
relationship between the level statistics and delocaliza- 
tion border has been much studied and remains in gen- 
eral an open question 0, H, 13 ■ 

In order to identify parameter ranges where chaos is 
present, it is necessary to quantify how accurately the 
statistical properties of the eigenvalues and eigenvectors 
of the Hamiltonian are described by RMT. Here, we 
will address this question by examining the energy level 
and eigenvector statistics across the full spectrum, rather 
than restricting to a specific spectral region (see also 
Sec. IV. C for additional quantitative discussion of this 
point). 

The extent to which the energy level spacing distri- 
bution distribution interpolates between the Poisson and 
Wigner-Dyson limits may be conveniently parametrized 
by a so-called Level Statistics Indicator (LSI), introduced 
in lii: 



/;°[Pp(s)-P„B(s)]ds 



(4) 



where so ~ 0.4729 is the first intersection point of Pf{s) 
and Pw>{s). LSI has the value 77 = 1 when P{s) = i-p(s), 
whereas = if P{s) = Pwd(s)- 

In order to characterize the degree of eigenvector delo- 
calization, a convenient measure is the so-called Number 
of Principal Components (NPC) [3, [13, [EH ■ Given a basis 
{\n)} of the system Hilbert space 7i, NPC of a normalized 
pure-state vector is defined as follows: 



(5) 



Thus, NPC estimates the number of basis states relative 
to which IV'} has a significant component. For example, if 
\ip) is a uniform superposition of exactly m basis states. 



then ^(|'0)) = '71. For an A'^-dimensional state vector 
with a component distribution pertaining to the GOE, 
the expected NPC is given by 



Cgoe 



where the factor 1/3 emerges from the Gaussian fluctu- 
ations of the eigenstates, and the additive correction 2 
is due to normahzation [41,152]. Since N is usually large, 
the approximation ^goe ~ N/2> is commonly adopted. In 
parallel with the fact that the LSI as constructed takes 
into account the level statistics across the full spectrum, 
we will be primarily interested in the average value of 
NPC across all of the eigenvectors in a relevant subspace. 



III. DISORDERED HEISENBERG MODELS 

The Heisenberg model in one- and two-spatial dimen- 
sions plays a paradigmatic role in condensed matter 
physics and statistical mechanics as a testbed for explor- 
ing quantum magnetism and spin dynamics in reduced 
dimensionality [53]. Thus, an accurate characterization 
of its integrability properties in physical regimes of in- 
terest has both a fundamental and practical significance. 
In the special case of spin-1/2 particles, a partial char- 
acterization of the integrability-to-chaos transition has 
been achieved, based on both clean systems where chaos 
is induced by coupling two different spin chains 1541 o r 
by adding next-nearest-ncighbor interactions (s^ . I55l.l56l] , 
as well as disordered systems, where the integrability- 
breaking term consists of random mag netic fields applied 
to all or a subset of spins [53, [H, H^j- Aside from their 
relevance to model real materials, disordered systems of- 
fer the added advantage of providing a natural arena 
to study the interplay between interaction and disorder, 
which remains a most challenging problem in condensed- 
matter physics. 

Within QIS, the isotropic Heisenberg spin-1/2 model in 
an external magnetic field arises naturally in some of the 
most promising solid-state proposals for scalable quan- 
tum computation, each spin corresponding, in the sim- 
plest settings, to a logical qubit and the exchange interac- 
tion providing the required inter-qubit coupling. Follow- 
ing the original suggestion by Loss and DiVincenzo [2^ 
for coupling electron-spin qubits via tunable Heisenberg 
interactions in semiconductor quantum dots, schemes for 
effecting exchange-based universal quantum computation 
have been further developed for both electron ^60j and 
donor-atom nuclear spins (30| . constructive methods for 
universal quantum gate design and efficient readout being 
identified in [61]. In addition, scalable universal archi- 
tectures where always- on Heisenberg couplings are used 
in conjunction with appropriate encodings of a logical 
qubits into three or more physical spins have been con- 
structed, offering both substantial implementation flexi- 
bility [29|] and enhanced decoherence suppression [s^, [12] ■ 
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Depending on implementation detail, imperfect qubit 
fabrication and/or uncontrolled residual spin-spin cou- 
plings during storage or gating may introduce an inte- 
grability breaking perturbation H' , causing the prereq- 
uisite mapping to well-defined logical qubits to be lost. 
While a number of error control schemes exist in principle 
to counteract the effects of H' (notably, dynamical refo- 
cusing methods for static disorder as considered here, see 
e.g. [63|), understanding the error behavior due to H' 
remains an important preliminary step. 

The Hamiltonian for an L-site lattice of spin-1/2 par- 
ticles coupled by the Heisenberg interaction and subject 
to a bias field in the z direction is given by: 



(6) 



where a^^^ is the vector of Pauli matrices {ux^ , Uy^ , cri*'' ) 
acting on the two-dimensional Hilbert space of the i-th 
site. The parameter determines the on-site Zeeman 
energy of the i-th spin, and will be parameterized as 
Ei = £ + 5£i, where e is a non-zero average and 5£i 
are uniformly distributed within [— (5e/2, fe/2], charac- 
terizing the different strengths of local random magnetic 
fields. The interaction strength between spin i and j 
is given by J^j, which will be likewise parameterized as 
Jij — J + SJij, where J is the average coupling, and 
dJij represents random interactions, being uniformly dis- 
tributed in [— (5J/2, SJ/2] [gl]. Since we are interested in 
the whole spectrum, the sign of the exchange coupling pa- 
rameter J is irrelevant. In the numerical simulations, we 
shall assume J > 0. The set of interacting pairs {i,j} 
is determined by the topology of the lattice. We shall 
consider nearest neighbor interactions on a ID chain and 
on a 2D rectangular lattice. In addition, the following 
notation is introduced: 



r(0 



1=1 



Hj^-Y, ^^^^ • ■ 



(7) 
(8) 



As discussed in Sec. IV, these two terms correspond to 
integrable limits of the Hamiltonian given in Eq. 



A. Two-body Random Ensemble 

It is important to appreciate that the class of systems 
considered here (Eq. ([6])) is more accurately described by 
the so-called two-body random ensemble (TBRE), which, 
instead of L-body couplings as implicit in the GOE, in- 
volves only two-body interactions. The TBRE, and more 
broadly the embedded Gaussian ensembles with fc-body 
interactions {k < L), was introduced in [gl, [6^ as a 
physically more realistic statistical setting for describing 



few-body interacting systems, such as atoms, molecules, 
and nuclei. Among the differences between the two en- 
sembles [63, m, [6^ , we highlight the ones which are most 
directly relevant to this work: 

• The GOE local density of states (also called the 
LDOS profile, see e.g. 0,[5|) as a function of en- 
ergy shows a semicircular law, whereas it is Gaus- 
sian for the TBRE. 

• The TBRE lacks ergodicity, in the sense that the 
statistical properties of each ensemble member need 
not coincide with the ensemble average [70j . 

• Contrary to the energy-independent NPC value pre- 
dicted by the GOE, the NPC estimate for the TBRE 
shows a strong variation with the energy of the 
state, approaching the GOE prediction of N/3 
mostly in the middle of the energy spectrum. 



B. Relevant Symmetries 

Although integrability may be regarded as equivalent 
to the presence of a c omp lete set of symmetries of the 



defining Hamiltonian [7l|, if such a set is incomplete 



then chaos may still be present within distinct symmetry 
sectors. However, a Wigner-Dyson distribution occur- 
ring independently in different sectors will tend to be 
washed out if the corresponding energy levels are mixed 
together. Thus, in order for clear conclusions to be drawn 
based on the energy level statistics, it is necessary to de- 
symmetrize the spectrum according to its trivial sym- 
metries, that is, the level spacing distribution must be 
examined separately in each resulting symmetry sector. 
The symmetries of the Hamiltonian in Eq. ([6]) are deter- 
mined by the choice of model parameters and Jij : 

1. Rotational symmetry around the z-axis. For all val- 
ues of £i and Jij, the z-component of the total 

spin, Sz = '^z'V^j is a conserved quantum num- 
ber, that is, [Sz,H] = 0. For even L, the largest 
subspace corresponds to Sz = 0, having dimension 
iVo = (^/J = ^!/[(-^/2)!]^- We shaU focus primar- 
ily on this subspace, TYq, as chaos sets in here first. 
A natural basis for Tig is the one associated with 
the eigenstates of Hz, which will be referred to as 
the computational basis or simply c-basis. 

2. Conservation of total spin. When the on-site en- 
ergies are degenerate, Si — e, the total spin, S"^ — 
(Si^^'V^)^) also becomes a good quantum num- 
ber, [S'^jH] = 0. Therefore, in this parameter 
regime the system must be separately studied in 
each symmetry sector characterized by fixed quan- 
tum numbers Sz and S [72|. 

3. Symmetries due to lattice geometry. For £i = e 
and Jij = J, symmetries under site permutations 
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related to the geometry of the lattice must also be 
considered. Under periodic boundary conditions, 
H is invariant under cyclic translations and reflec- 
tions of the lattice sites, resulting in momentum 
and parity conservation respectively. Note that in 
rectangular 2D lattices, momentum and parity in 
both directions are good quantum numbers. For 
open boundary conditions, the Hamiltonian is in- 
variant under lattice reflections only. We shall re- 
strict our analysis to open boundary conditions in 
what follows, because they lead to fewer symme- 
tries hence larger invariant subspaces, which per- 
mit better eigenvalue statistics. Therefore, when 
identifying invariant subspaces, the eigenstates will 
be grouped according to their quantum numbers 
{Sz ,S,R), where R indicates parity. 



4. Time-reversal invariance. GOE is the ensemble ap- 
propriate for describing the static properties of sys- 
tems exhibiting time-reversal invariance, that is, for 
Hamiltonians that commute with the time-reversal 
operator, Tq = (i)^ (E)f=icry^ K , where K is the con- 
jugation operator in the c-basis for spin-1/2. It is 
important to note that the Heisenberg model with 
a magnetic field does not commute with the conven- 
tional time reversal operator, Tq. Nevertheless, the 
GOE rather that the more general Gaussian uni- 
tary ensemble (GUE) is appropriate for this model 
because the Hamiltonian matrix admits a straight- 
forward real representation in the c-basis. Equiva- 
lently, one may understand the applicability of the 
GOE as resulting from the fact that the Hamilto- 
nian commutes with K which may be regarded a 
non-standard time-reversal operator as in ^41. . .57i] . 



C. Related Models 

The Heisenberg model may be regarded as one of a 
related class of spin-1/2 lattice models with two-body 
interaction Hamiltonians of the form: 

Hmz = XI ^^'^^ 

i 

In particular, the 2D Ising model in a transverse field 
(Jf. = Jf. = 0), as well as the XXZ model (Jf^- = 
Jfj, Jfj ^ 0) and the isotropic XY model {J[j = Jf^, J^j = 
0), have been previously studied in the context of quan- 
tum computation and quantum chaos 0, 
[76| . However, the Heisenberg model differs in some es- 
sential characteristics, notably the nontrivial role of sym- 
metries. Also distinctive of the Heisenberg exchange cou- 
pling is the competing nature of the two interacting terms 
in the Hamiltonian: in the c-basis, the diagonal Ising in- 
teraction ai*'' favors localization, whereas the flip-flop 
term cri Vi + <^y "y induces delocalization. Thus, a 
comparison between our results and those of previous 
studies shall also serve to identify those properties which 
are generic to disordered spin-1/2 lattices as contrasted 
to those specific to the Heisenberg model. 

IV. RESULTS: LEVEL STATISTICS AND 
DELOCALIZATION PROPERTIES 

We shall proceed by first obtaining a basic characteri- 
zation of the level spacing distribution and average delo- 
calization properties of the eigenvectors for a wide range 
of parameters of interest, and then proceed to examining 
entanglement properties in Sec. VI. 



Although, as mentioned, for a fully established (or 
weakly broken) symmetry one may examine the level 
statistics in each invariant subspace, for strongly bro- 
ken symmetry it is no longer possible to sort eigenvalues 
by symmetry quantum numbers nor can correlations be- 
tween different subspaces be ignored. As a result, sepa- 
rating the effect of symmetries on the energy level spac- 
ing distribution from that of integrability vs. chaotic- 
ity becomes challenging throughout the symmetry tran- 
sition regime, that is, it is difficult to distinguish sys- 
tematically between the emergence of a complete vs. in- 
complete set of symmetries. While some analytical re- 
sults exist for eigenvalue statistics across the full space 
in the presence of symmetry-breaking ^73|, JA] , they do 
not offer a practical diagnostic tool for distinguishing 
between fully chaotic and partly integrable behavior in 
broken-symmetry regimes based on LSI (especially when 
multiple symmetry subspaces are involved as in the S'^- 
symmetry breaking which plays a relevant role here). 



A. From Integrability to Chaos 

The model described by Eq. ^ shows two limiting 
integrable cases, Hz and Hj [Eqs. ([7]) and ([8]), respec- 
tively]. The first is a trivially solvable non- interacting 
problem, while the second is solvable only in ID by Bethe 
Ansatz [tJ [zl, [7§, Hoj . A transition to chaos may be in- 
duced by adding an integrability-breaking term to any 
of the two integrable limits. Here, we shall focus on the 
following representative scenarios: 

Case J/Se: Exchange interactions with a constant 
strength J are added to Hq = Hz- The crossover 
between integrability and chaos is studied as a func- 
tion of the ratio J/Se, with 5J = throughout. 

Case (5 J/ J; Exchange interactions with random 
strength are added to Hq = Hj. The crossover 
between integrability and chaos is analyzed as a 
function of the ratio SJ/J, with & = throughout. 
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Case 5J/Se: Exchange interactions with random 
strength are added to Hq = Hz- The crossover 
between integrabiUty and chaos is studied as a 
function of the ration SJ/Se. 

Clearly, Case J/6e includes the other two cases in 
limiting situations: Specifically, it coincides with Case 
6 J I Se when the on-site disorder dominates over the dis- 
order in the coupling strengths, 6e — > oo, and with Case 
6 J I J when the distribution of coupling strengths is very 
narrow, J — > oo. 

B. Delocalization Measures 

Because the properties of a state vector (in particular, 
delocalization) depend entirely on the choice of represen- 
tation, quantities such as NPC cannot serve as intrinsic 
(basis-independent) indicators of quantum chaos, such 
as energy level statistics are considered to be. This fact 
may undermine the entire program of establishing con- 
nections between the properties of the eigenvectors of a 
Hamiltonian and its level of chaoticity. Yet, it is well- 
known that for chaotic Hamiltonians, RMT predictions 
hold (to some extent) in a large number of reasonable 
representations, although a systematic characterization 
of such representations and the degree of agreement one 
may expect remain at present an unanswered question. 
Here, we conform to the standard approach to this prob- 
lem, and examine eigenvector properties with respect to 
a basis in which the spread of the eigenvectors has a phys- 
ical meaning, for instance one which relates to the rele- 
vant measurement capabilities, or to an integrable limit 
of the class of Hamiltonians under consideration. On one 
hand, we shall investigate to what extent the eigenvector 
properties associated with quantum chaos depend on the 
choice of two different bases corresponding to the inte- 
grable limits (IZl)-®- On the other hand, we shall intro- 
duce a measure of relative delocalization between bases 
associated with different disorder realizations. Specifi- 
cally, we will examine: 

(i) The above-mentioned c-basis^ which is associated 
with the eigenstates of Hz- In quantum compu- 
tation, this is the basis that represents classical in- 
formation, relative to which the final read-out mea- 
surement is performed The associated NPC will 
be labeled £,c- 

(ii) A so-called J-basis, or interaction basis, which cor- 
responds to the eigenbasis of Hj. In this case, the 
associated NPC will be denoted by ^j. 

(iii) A disorder-dependent relative representation, which 
is obtained as follows. Within a set of sequentially 
generated random realizations, the eigenbasis of a 
given random realization is used to calculate NPC for 
the eigenvectors of the subsequent realization. The 
associated NPC, which will be denoted ^r, quantifies 



relative delocalization^ that is how delocalized are 
the eigenbases associated with different disorder re- 
alizations with respect to each other, rather than 
with respect to some fixed basis. By construction, 
this quantity is independent of any fixed basis and 
depends only on the disorder. 

Interestingly, an approach which attempts to quantify 
the "complexity" of the eigenstates of an ensemble of 
Hamiltonians in a way similar to the above-mentioned 
relative delocalization has been jiroposed based on the 
notion of correlational entropy [81]. Both correlational 
entropy and relative delocalization overcome the prob- 
lem of basis dependence by invoking a distribution over 
Hamiltonians. An important dilference, however, is that 
evaluating the correlational entropy requires tracking in- 
dividual eigenstates as a function of the disorder across 
the full disorder space. While a further comparison be- 
tween the two concepts would be worth pursuing, it is 
likely that this feature could make relative delocaliza- 
tion computationally more tractable for many disorder 
settings. 



C. Numerical Results 

We consider the ^2 = 0-subspace, Ho, for ID and 
2D models. In both cases, for lattice size L = 12, 
dim (TLo) = Nq = 924. As anticipated in Sec. II, NPC val- 
ues are averaged over all eigenstates in Ho, and both LSI 
and NPC are further averaged over a number of disorder 
realizations sufficient for accurate statistics, as indicated 
in each figure caption. 

Case J/6e. For this disorder setting, two 
integrability-chaos transitions are verified for the ID 
model, whereas only one occurs in 2D. From the top pan- 
els of Fig.[l] we see that in both cases the first crossover, 
from the Hz integrable limit to chaos, is observed as 
the interaction strength J increases from zero to a value 
close to the energy difference between directly coupled 
states, J/Se ~ 1. In the other extreme, where J/fe — > oo 
and therefore H Hj, the transition to 5*^- and R- 
symmetry complicates the interpretation of the energy 
level statistics. During this transition, different {S, R) 
invariant subspaces partially overlap, and the observed 
LSI increase across the whole Sz = 0-subspace need not 
reflect a change towards integrability, but rather the pro- 
gressive decoupling of states belonging to difi^erent sub- 
spaces (this region is indicated by a dashed line in the 
figure). Once the transition is complete, at J/Se —^ oo, 
the level spacing distribution within individual {S^ , S, R) 
symmetry sectors may be determined. This results in 
a Poisson distribution for the ID model [s^l, consistent 
with its exact solvability. However, and contrary to the 
behavior of the LSI across the whole Ho, the level spacing 
distribution in the invariant subspaces of the 2D system 
are strongly Wigner-Dyson. 

It is interesting to contrast the LSI results with those 




FIG. 1: (Color online) Case J /5e. 5*2 = subspace, L — 12. 
Top panels: Level statistics indicator, rj, vs. J/5e. The 
dashed line indicates in each case the transition region, where 
the {S, R) subspaces are not fully distinct. The arrow indi- 
cates the actual limiting value of the LSI value as J/5e ^ oo. 
Bottom panels: Number of principal components, ^, vs. J/Se. 
Left panels: ID chain, right panels: 2D 3x4 lattice. Averages 
over 20 random realizations. 



FIG. 2: (Color online) Case SJ/J. (S^,S) = (0,1) sub- 
space, L = 12. Top panels: Level statistics indicator, rj, vs. 
SJ/J. Dashed lines indicate in each case the region where the 
transition to J?-symmetry sectors occurs. Dotted line: Gaus- 
sian disorder with a = 5 J /A. The arrow indicates the actual 
limiting value of as SJ/J —> 0. Bottom panels: Number of 
principal components, vs. SJ/J. Left panels: ID chain, 
right panels: 3x4 lattice. Averages over 50 random realiza- 
tions. 



for NPC. The dependence of £_c, ^-nd £,r on J/6e is il- 
lustrated in the bottom panels of Fig. [TJ Corresponding 
to the broader chaotic region detected by the LSI, there 
is greater delocalization in 2D, especially as quantified 
by and ^r-. For both cases, the maximum NPC in most 
instances occurs within the chaotic region, which we take 
to be ?7 < 0.3 [7]. However, the actual value falls short 
of the RMT prediction and different NPC measures may 
disagree substantially even within the chaotic region. For 
example, at J/Se « 2 where reaches its maximum, 
and £,r show only partial delocalization. and de- 
tect only one localized-to-delocalized transition, whereas 
is only large where LSI is low. Note that the behavior 
of matches more closely the behavior of LSI on the 
low J/Se-LSl transition than ^c- However, in the limit 
J/de —> oo fails to discriminate between the integrabil- 
ity of the ID model and the presence of chaos in the 2D 
model, decreasing to 1 in both cases due to the absence 
of disorder. 

Case SJ/J. In this case, [S'^,H] — 0, thus we exam- 
ine LSI and NPC in the largest S'^-subspace, correspond- 
ing to 5 = 1 and dimension N = 297 eigenstates. As 
6 J/ J {H — > Hj), the reflection symmetry becomes 
important, causing the LSI in the top panels of Fig.[2]to 
increase irrespective of the chaoticity or integrability of 
the Hamiltonian. As previously discussed, in this limit 
the ID model is integrable, while the 2D model is not. 
Consequently, as dJ/ J increases from 0, a transition from 



integrability to chaos occurs for the ID model, whereas 
only the breaking of the reflection symmetries occurs in 
the 2D model. Notably, for the ID model near 6 J/ J — 2, 
an unexpected rapid transition in LSI from r/ ^ 0.2 to an 
intermediate value of 77 ~ 0.6 is observed. Although the 
2D system remains chaotic throughout, a small rise in 
LSI close to 5 J/ J ^ 2 is noticeable. 

NPC values are shown in the bottom panels of Fig. [21 
Since the Hamiltonians in the J-basis and in any basis 
associated to a disorder reahzation are block-diagonal, 
and cannot exceed the dimension of the {Sz, S) — (0, 1) 
subspace, unlike which is upper-bounded only by the 
dimension of the Sz = subspace given that states of the 
c-basis do not posses S'^-symmetry. For the ID model, 
the slight increase in as 6 J/ J 1 occurs during the 
integrability-chaos transition, but the same occurs also 
for the 2D system, which is chaotic throughout. Thus, 
this increase is likely related to the reflection symmetry- 
breaking. Interestingly, the abrupt drop in for both 
ID and 2D near SJ/J = 2 appears to be connected with 
the rise in LSI - the greater change in both LSI and 
occurring for the ID model. 

A possible explanation for the observed transition at 
6J/J ~ 2 is based on the fact that this value marks 
the point where interaction strengths arbitrarily close to 
zero are allowed, mmij{Jij} — J — 5 J/2 = 0. If the 
chain is broken into two or more approximately uncou- 
pled segments, then the energy spectrum of each segment 
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becomes approximately independent. Even if for each 
segment the level spacing distribution is Wigner-Dyson, 
the level spacing distribution for the combined spectrum 
will not be. Lending support to this explanation is the 
far stronger effect in ID than in 2D, following from the 
fact that the 2D lattice cannot be as easily broken into 
isolated segments. In order to determine if this effect is 
an artifact of the sharp cut-off in the disorder distribu- 
tion, interactions with Gaussian disorder (with standard 
deviation a = SJ/4) were also examined. The results 
are shown as the dotted line in Fig. [2] (top left), where 
it is seen that the abrupt transition remains despite the 
absence of a sharp cut-off in the disorder range (note, 
however, that the width of the transition is affected). 

It is intriguing to notice that in ID the ground state 
and low-energy excitations for this disorder setting have 
been determined at 6 J/ J = 2 through an exact renormal- 
ization group approach, and are known to be described 
by the so-called random singlet phase [ssl . Is^ [ssj . A 
more thoroughly discussion of possible connections be- 
tween the observed behavior and the existence of a ran- 
dom singlet phase are left for future investigation. 

As a further remark, the behavior of in the region 
where SJ/ J > 2 is worth mentioning. In ID, S.r is much 
lower than ^j, while in 2D only a small difference is seen. 
This further demonstrates that, in the presence of dis- 
order, relative delocalization can be an effective indicator 
of chaos. 

Case dJ/Se. The dependence of LSI on 6J/5e is 
shown in the top panels of Fig. [3] As the disorder in the 
coupling strengths becomes sufficiently large to compete 
with the energy difference Se, a transition from integra- 
bility to fully developed chaos is observed only for the 
2D model, whereas in ID theLSI merely approaches an 
intermediate value of 77 ^ 0.6 ^8d\. For SJ/5e > 10, LSI 
for the entire Sz — subspace increases, but, as before, 
this does not necessarily reflect an approach to integra- 
bility, since a transition to the S^-symmetry is in place. 
In fact, at SJ/Se — > 00, we have verified that essentially 
the same LSI values reached at SJ/6e ~ 10 for ID and 
2D within the {Sz,S) = (0, 1) subspace are attained. 

The NPC behavior is depicted in the bottom panels of 
Fig. [3l Refiecting the fact that chaos never fully devel- 
ops in ID, both and S.r are considerably lower when 
compared to the previous case. In contrast, for the 2D 
model, this disorder setting leads to a high level of delo- 
calization, especially in the J-basis and in terms of rel- 
ative delocalization. The maximum values are reached 
between 1 < SJ/Se < 6. As SJ/Se increases further, the 
transition to S'^-symmetry sectors becomes relevant, and 
the Hamiltonians in these two bases approach a block- 
diagonal form - therefore and decrease accordingly. 
As SJ/Se 00, the values obtained still indicate strong 
delocalization relative to the sizes of the (0, 5") subspaces. 

Consistent with studies of the TBRE and of other sys- 
tems with two-body interactions (notably, nuclear-shell 
models), a strong dependence of NPC on energy is ob- 
served, particulary in chaotic regions, see Fig. 31 Specif- 
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FIG. 3: (Color online) Case SJ/Se. Sz = subspace, 
L = 12. Top panels: Level statistics indicator, ri, vs. SJ/Se. 
Dashed line indicates the region where the transition to S^- 
symmetry sectors occurs. Arrow indicates the actual limiting 
value of ?7 as SJ/Se 00. Bottom panels: Number of prin- 
cipal components, ^, vs. SJ/Se. Left panels: ID chain, right 
panels: 3x4 lattice. Averages over 50 random realizations. 



ically, the NPC follows an approximately Gaussian form. 
We show a representative 2D case at SJ/Se ^ 3, where 
f J and are close to their maximum. Note that RMT 
provides a bound on the extent of delocalization, and 
that even in the center of the spectrum, the average de- 
localization falls short of the predictions of RMT, partic- 
ularly in the c-basis. That the eigenstates in the tails of 
the spectrum are less delocalized than those in the center, 
however, only partially affects the average NPC's shown in 
Figs. [11 [21 [21 because the density of states is also peaked 
around the center of the spectrum. We found that aver- 
aging only over the central 100 eigenstates (10.6% of the 
total) increases the values of each NPC -measure by an 
amount between 10% to 15%, and does not qualitatively 
affect the behavior of the average NPC. For all cases ex- 
cept for and S.r in case SJ/Se near their maxima, the 
deviation from the RMT prediction of remains sub- 
stantial. 

This situation is to be compared with the models stud- 
ied in ^ and '26*1 . Even though these models also corre- 
spond to the TBRE as well, better agreement with RMT 
is obtained for delocalization whenever the interactions 
are either purely off-diagonal or all-to-all, that is, ev- 
ery site is coupled with every other site. In fact, a key 
difference in our case is the presence of the purely diag- 
onal Ising contribution, which favors localization in the 
computational basis. Furthermore, we have verified that 
all-to-all Heisenberg couplings also yield delocalization 
values in better agreement with RMT (data not shown) . 
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FIG. 4: (Color online) Case dJ/Se. Number of principal 
components vs. energy eigenvalue for SJ/5e = 3. Eigenstates 
in the Sz ~ 0-subspace for a 3 x 4 qubit lattice are considered. 
The dashed horizontal line shows the GOE prediction, ^goe = 
(A^o + 2)/3 308. Averages over 20 random realizations. 

Thus, the lower level of ergodicity in the models exam- 
ined here may be partly attributed to lower connectivity. 



V. ENTANGLEMENT MEASURES 

Obtaining a complete characterization of entanglement 
in many-body systems is a challenging problem which in 
spite of extensive effort remains as yet largely unsolved 
hj. As a main reason for such difficulty, it is evident 
that no single entanglement measure can fully capture 
the complexity of multi-particle correlations. For our 
current purposes, we shall select representative measures 
of entanglement that are both computationally practi- 
cal, and directly connect with prior work on entangle- 
ment and quantum chaos and/or quantum phase tran- 
sitions. Specifically, we shall investigate so-called con- 
currence between selected pairs of qubits as a measure 
of pairwise correlations [87|, and a family of multipar- 
tite purity measures constructed within the general GE 
framework [36l.[37j. 

For a pure state lip) of two qubits, concurrence may be 
defined as 

where 1-0*) is the complex conjugate of Physically, 
C{\ip)) may be thought of in this case as the overlap 
of a state with its time-reversed counterpart. When a 
pure state lip) of L > 2 qubits is considered, the reduced 



two-qubit state of a selected pair {i,j) is described by a 
density matrix pij , where all but the qubits of interest are 
traced out. In this case, concurrence may be computed 
through a more general expression, which holds also for 
mixed states and reads 

) = max(0, Ai - A2 - A3 - A4), (9) 

where Ai > A2 ^ A3 ^ A4 are the square roots 
of the eigenvalues of PijPij, and pij = (cy^'' O a-y^^) 

Piji^v^'' ® CTy^^) [13 ■ Concurrence ranges from a mini- 
mum of for unentangled states to a maximum of 1 for 
states containing a maximum amount of pairwise cor- 
relations. Thus, applied to a pair of qubits in a pure 
many-qubit state, zero concurrence will occur for both a 
product state, lip) = (8) ... €5 but also for the 

L-partite Greenberger-Horne-Zeilinger state, |GHZi) — 
-^(|0,...,0) -I- which exhibits genuinely mul- 

tipartite correlations For random states of a suf- 

ficiently large number of qubits, the expected value of 
concurrence between any two qubits is zero [s^. 

For a multipartite system, concurrence quantifies the 
amount of mixed-state bipartite entanglement within a 
given pair. Complementary information about how the 
pair itself or, more generally, a given subset A of qubits, 
correlates with the remaining subset i?, is provided by the 
amount of bipartite entanglement between A and B. The 
unique measure satisfying all requirements of invariance 
under local transformations, continuity, and additivity 
is the von Neumann entropy of either reduced density 
matrix, e.g., 

S{\iP)ab) = -Tta{pa \og2 pa), 

where pA is the reduced density operator of subsystem 
A. If the additivity requirement is relaxed, a simpler lin- 
earized version of the above expression suffices to quan- 
tify bipartite entanglement, the so-called linear entropy, 

E{\^)ab) = 1-Tta{pa)^. 

Although the amount of multi-partite entanglement is 
generally hard to quantify away from bipartite states, 
useful indirect insight may be gained by considering dif- 
ferent bi-partitions of the system. In particular, \ip) cer- 
tainly contains genuine multipartite entanglement if no 
reduced subsystem state is pure. 

The notion of GE offers a powerful framework for 
both organizing conventional multipartite entanglement 
within a unified setting, and for extending the concept 
of entanglement to situations where a preferred partition 
of the system into subsystems may not be meaningful or 
otherwise desirable [3y, [33, [3^, [90] ■ GE is based on the 
relationship of the state of interest to a distinguished set 
of observables, rather than to a tensor product decompo- 
sition of the Hilbert space Ti into subsystems. Let such 
a distinguished observable set consists of the Hermitian 
operators in a linear subspace h of the full operator space 
on 7i, with h closed under Hermitian conjugation. The 
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key step is to replace the notion of subsystem state as ob- 
tained via the usual partial trace operation by a notion of 
"reduced state" as determined by the expectation values 
of observables belonging to the restricted subspace h. A 
pure state \ip) is generalized unentangled (entangled) rel- 
ative to h depending on whether its reduced state is pure 
(mixed) in the space of all reduced states that is, depend- 
ing on whether it is extremal (or not) in the convex sense. 
If {Bi] is a basis of Hermitian traceless operators for ft,, 
orthogonal in the trace inner product, a natural measure 
for GE is the degree of purity of the reduced /i-state as 
quantified by the so-called h-purity, 

i 

where the overall normalization constant k is chosen so 
that Ph ranges between its maximum value of 1 for gen- 
eralized unentangled states and its minimum value of 
for maximal GE relative to h. 

When the observable set h is the Lie algebra of all lo- 
cal observables on qubits, GE reduces to standard mul- 
tipartite qubit entanglement. In particular, the above h- 
purity (fTU)) relative to arbitrary single-qubit observables 
takes the explicit form: 

i=l,L 

Piim^Y ^ KV'ki^^iV')!'. (11) 

a—x.y,z 

This quantity is simply related to the so-called Meyer- 
Wallach measure of global entanglement, Q, by Pi = 1 — 
Q 38, 91, 92], which quantifies multipartite entanglement 
through the average bipartite entanglement between a 
qubit and the rest. By construction, Pi = 1 for product 
states, whereas Pi — for states such as |GHZi), where 
each single-qubit reduced density matrix is totally mixed. 

A family of related entanglement measures may be nat- 
urally constructed from the above local purity by coarse- 
graining. If the qubits are partitioned into distinct n- 
qubit blocks, the n-local purity, denoted Pn, is defined as 




FIG. 5: (Color online) Choices of partitions for evaluating 
n-local purities in the 3x4 lattice. 



Pi P2 Pa Pa Pe 
2.16 5.69 7.71 7.49 10.10 

TABLE I: Expected n-local purity for random states in the 
Sz = Q subspace of a L = 12-site lattice, A^'o = 924. Values 
are multiplied by a factor x 10^ for clarity. 

purity with respect to arbitrary observables local to each 
block. For example, the bi-local purity (also used in psj ) 
is given by 

^2(1^)) = ^ 'e' m-^^'-'~'-?m^ (12) 

a.j3—x.y,z 

where the sum only extends to traceless operators. Phys- 
ically, the measure Pn, for n > 1, ignores short-range 
correlations which Pi detects. For the ID chain we con- 
sider, the natural choice of partitions is into contiguous 
blocks. For the 2D lattice, relevant block partitions are 
illustrated in Fig. [5] 

It is important to mention that the von Neumann en- 
tropy of an n-qubit block, which has been a more broadly 
used quantity in studies of entanglement vs quantum 
chaos 19, 20, 26, 27], is closely related in meaning to P„. 
It can be shown that the average bipartite entanglement 
between the n-qubit block and the rest, Q„ = 1 — P„, is 
(up to a normalization constant) equal to the average lin- 
ear entropy over each n-qubit block participating in P„. 
Thanks to the average over many blocks, P„ is less sensi- 
tive to edge effects. Apart from that, a main advantage 
with respect to von Neumann block entropy is its math- 
ematical simplicity, which allows analytic calculations of 
the expected GE for random pure states to be and also 
reveals a quantitative relationship with delocalization as 
measured by NPC in a local basis jl^]. In particular, for 
random pure states with purely real components, as ex- 
pected for GOE eigenvector statistics, the average value 
P„ with respect to the appropriate Haar measure may 
be exactly computed for arbitrary lattice size L. While 
we defer the details of the derivation to Appendix A, the 
results for states in the Sz = 0-subspace with L = 12 are 
summarized in Table I. 

The intrinsic flexibility of the GE notion easily allows 
for the construction of purity measures relative to distin- 
guished operator spaces not directly tied to a subsystem 
partition. For spin-1/2 systems, in particular, a natural 
choice of distinguished observables emerges upon map- 
ping Pauli spin operators to canonical fermionic opera- 
tors via the Jordan- Wigner transform. 

It is straightforward to show that the {ci} satisfy canon- 
ical fermionic anti-commutation relations, {cJ,Cj} = 5ij, 
{ci,Cj} — 0. Then "generalized local" resources may 
be associated with quadratic fermionic observables com- 
muting with the total fermionic number operator, L = 
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cjci - as opposed to "nonlocal" resources involved in 
processes where the total fermion number in a given pure 
state may change. The corresponding distinguished ob- 
servable set is isomorphic to the unitary Lie algebra in 
L dimensions, h = u{L). GE relative to such algebra is 
quantified by the fermionic u(L)-purity [s^, 



L 



(13) 



Physically, Pu{l) indicates how "close" a state is to be- 
ing described by a fermionic product state (a Slater de- 
terminant) ^9^. For example, in the two excitation sec- 
tor, Pu(L){\4')) = 1 if and only if \ip) may be written 
in the form ~ c|jcj^|vac), where a,b label any set of 
modes unitarily related to modes and |vac) contains 
no fermions. Pu(l) has been shown to successfully de- 
tect and characterize broken-symmetry quantum phase 
transitions (ssj . In the context of the transition to quan- 
tum chaos, such a measure may be expected to provide 
insight into entanglement generation associated with the 
departure from a non-interacting fermion problem. 



VI. RESULTS: ENTANGLEMENT BEHAVIOR 

We are now in a position to address entanglement prop- 
erties of Heisenberg-model eigenvectors, and to under- 
stand the resulting behavior based on the insight gained 
from the corresponding behavior of delocalization and 
energy level statistics. 



A. Disorder Dependence 

Case J/6e. As remarked, this setting is espe- 
cially useful for exploring possible connections between 
integrability-breaking, delocalization, and entanglement, 
because it contains both the integrable limits associated 
with Hz (in ID and 2D) and Hj (in ID). With respect 
to the c-basis, the ID model possesses three regimes of 
interest: (i) integrability-localization for J/Se 0; (ii) 
chaoticity for intermediate ratios J/Se ^ 1; and (iii) 
integrability- delocalization for J/Se ^ oo in ID, while 
in latter regime the 2D model remains strongly non- 
integrable. 

Fig. [6] summarizes the behavior of the average con- 
currence between all neighboring spins, C, and of the 
n-local purities, P„, for the ID and 2D models. As we 
depart from the purely local Hamiltonian Hz and SJ/ 6e 
increases from to 1, a peak in C is observed before the 
onset of chaos. This suggests that, as the system delocal- 
izes, pairwise entanglement is the first type of entangle- 
ment to emerge - disappearing, however, in the chaotic 



region [25|. Such a disconnection between pairwise en- 
tanglement and the onset of chaos has also been verified 
for the 2D Ising model in a transverse field in [2^ and, 
more recently, 27[. Contrary to that, each of the P„ de- 
creases from nearly 1 to nearly 0, roughly corresponding 
to the transition in LSI and (cf. Fig. [T]). Therefore, 
as we approach chaoticity, a shift from pairwise to gen- 
uinely multipartite entanglement occurs. Interestingly, 
in ID a shift from short-range to long-range correlations 
is also noticeable, since a more pronounced decay is wit- 
nessed by Pi , which is subsequently followed by the other 
purities until finally being reached by Pq. In 2D, differ- 
ent Pn do not directly signify correlations over different 
distance, and all P„ curves superimpose more closely. 

In the region where J/Se increases from 1 to oo and 
the limiting Hamiltonian Hj is approached, C and Pi 
are near zero throughout - note that due to the rota- 
tional invariance of Hj, Pi — > as J/Se — > c», whereas 
P2,P3,P4, and Pg increase slightly, particularly in ID. 
It is hard to ascertain whether this subtle variation in 
ID is related to a transition to integrability or it is sim- 
ply associated to the symmetry transition. Therefore, 
the similar entanglement behavior observed in both ID 
and 2D does not mirror the differences in chaoticity be- 
tween the two systems. In fact, it appears much more 
closely connected to the behavior in this region. In 
this sense, a parallel between standard multipartite en- 
tanglement and delocalization is stronger than between 
the former and chaos. Furthermore, similar to delocal- 
ization, the minimum values reached by each P„ in the 
central chaotic region, although low, do not attain the 
values predicted by RMT. For the 2D case, where agree- 
ment with RMT is best, the minimum of each P„ except 
for Pi range between a factor of 4.5 to 5.8 times the 
RMT values in Table I. Pi does attain the RMT value 
but only well into the S'^ symmetry region, where this 
effect cannot be attributed to chaos. 

Case dJ/J. Similar to the above Case J/Se in the 
region J/Se > 1, the entanglement behavior shows no 
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FIG. 6: (Color online) Case J/Se. Average nearest-neighbor 
concurrence and n-local purity averaged over all eigenvectors 
of the Sz ~ 0-subspace vs SJ/5e. Lattice size L = 12. Solid 
black line: concurrence. Left panel: ID chain, right panel: 
3x4 lattice. Averages over 20 random realizations. 
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qualitative difference between tfie ID and 2D models for 
5J/J G [0,1], although an integrability-chaos transition 
occurs in ID (see Fig. [71). On the other hand, both _P„ 
and C change abruptly near 5J / J 2 in ID, and near 
5J/ J ~ 3 in 2D. In both cases, the sharp increase in 
pairwise entanglement and decrease in multipartite en- 
tanglement are paralleled by a decrease in and, more 
interestingly, by an increase in LSI (cf. Fig. [2]). It is ex- 
pected that the breaking of the chain into distinct sub- 
systems will result in a decrease of multipartite entangle- 
ment. It is, however, intriguing that although smaller, a 
decrease occurs also in 2D, where chain breaking should 
be suppressed. 

Case 5J/6e. The entanglement behavior in 2D for 
this disorder setting is very similar to the one depicted 
in the right panel of Fig. [51 whereas in ID purities never 
reach low values, reflecting the lack of delocalization (cf. 
Fig. 121) . We shall then focus on the 2D system and de- 
scribe some interesting features of purity, which also ap- 
ply to Case J/Se in both ID and 2D. 



B. Energy and Delocalization Dependence 

In addition to understanding how entanglement prop- 
erties depend on disorder and interaction strength, their 
direct dependence upon energy and delocalization is rel- 
evant to gaining a more complete physical picture. 

Similarly to the behavior observed for NPC within the 
TBRE (see Fig. 01), a strong correlation also exists be- 
tween multipartite entanglement and the energy spec- 
trum. The left panel of Fig. [51 contains a representative 
example for the bi-local purity. In general, states at the 
edges of the spectrum tend to be less entangled, whereas 
highly entangled states are clustered around intermediate 
energies. This result indicates a correspondence between 
Pn and the LDOS of the TBRE, which, as a function of 
energy, is Gaussian and broadly peaked at the center of 
the spectrum [g^, [H, [6§| (see also HH for a discussion 
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FIG. 7: (Color online) Case 6 J/ J. Average nearest-neighbor 
concurrence and n-local purity averaged over all eigenvectors 
of the (Sz, S) — (0, 1) subspace vs 5 J/ J. Lattice size L = 12. 
Solid black line: concurrence. Left panel: ID chain, right 
panel: 3x4 lattice. Averages over 20 random realizations. 
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FIG. 8: (Color online) Case SJ/Se. Sz ~ subspace, 3x4 
2D spin lattice, SJ/Se = 2. Left panel: P2 vs. energy eigen- 
value. Right panel: Pi vs ^c- A single disorder realization is 
considered. 



of the influence of LDOS properties on purity behavior). 
Essentially, where the LDOS is largest, a higher level 
of delocalization exists and also higher amounts of mul- 
tipartite entanglement occur. We emphasize, however, 
that this relationship is not per se indicative of chaos: A 
similar behavior was found for disordered on-site energies 
(Case J/Se) within the {Sz,S,R) — (0,1, +1) sector for 
the ID model in the limit J/Se — + 00, which corresponds 
to a delocalized, but integrable regime. 

The relationship between multipartite entanglement 
and delocalization may further be probed by directly 
comparing n-local purities and NPC. A plot of the local 
purity. Pi, as a function of delocalization in the c-basis 
(right panel of Fig.[8l) discloses a striking relationship be- 
tween Pi and ^cj which is found to persist for each one 
the Pn and for a broad range of values of SJ/ Se within the 
chaotic region (until the S'^-symmetry becomes strong). 
Remarkably, for a given P„ , the precise shape of the curve 
does not depend on SJ/Se, although it depends slightly 
on lattice dimension and disorder setting. The precise 
relationship between Pi and is examined at length 
in [93| . where it is shown that the relationship between 
Pi and depends strongly on how correlated the state 
vector components are with respect to the Hamming dis- 
tance between the quantum numbers describing c-basis 
states 9J|. If no such Hamming-correlation exists, in- 



verse proportionality between local purity and NPC is ex- 
pected (see also 23]): 
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with = A^o = 924 for the central band. While this 
qualitatively agrees with the observed behavior, due to 
the fact that all interaction terms in the Hamiltonian are 
of a two-body nature, the components of the eigenvectors 
tend to be Hamming-correlated, resulting in significant 
deviations from the predicted inverse scaling law, espe- 
cially at small ^c- This effect was recently independently 
confirmed by Giraud and coworkers [951] . 

Beside calling for a deeper understanding of the phys- 
ical conditions leading to the observed non-trivial eigen- 
vector structure, the above findings naturally prompt the 
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FIG. 9: (Color online) Multipartite entanglement vs. delo- 
calization in a clean Heisenberg Hamiltonian: Pg vs in 
the {Sz,S,R) — (0, 1, +1) symmetry subspace. Lattice size: 
L = 12. Left panel: ID chain. Right panel: 3x4 2D lattice. 



FIG. 10: (Color online) Generalized fermionic entanglement 
vs delocalization in a clean Heisenberg Hamiltonian: Pu(l) 
vs. in the {Sz,S,R) = (0, 1, +1) subspace. Lattice size: 
L = 12. Left panel: ID chain. Right panel: 3x4 2D lattice. 



following question: To what extent could the relationship 
between entanglement and delocalization provide a signa- 
ture of quantum chaos? As a first step toward answering 
this question, we revisit the case of a clean Heisenberg 
Hamiltonian, Hj (Case J/5e in the limit of no disorder), 
which supports both integrability (in ID) and chaoticity 
(in 2D). Instead of Pi, which is identically zero for eigen- 
vectors oi H = Hj, we look at the block-purity Pg- As 
illustrated in Fig. [51 no clear relationship between and 
P„ emerges in the integrable regime (left panel) , whereas 
a noticeable relationship is present in the chaotic case 
(right panel) - in spite of the limited statistics accessi- 
ble due to symmetry constraints. Interestingly, this is 
in agreement with a recent study of entanglement and 
chaos for the Ising model in a transverse field by Lak- 
shminarayan and coworkers [2^ , where a correlation be- 
tween values of delocalization and block von Neumann 
entropy is found to occur only in chaotic regimes and 
not in integrable yet delocalized ones. 

While the above results provide suggestive evidence 
in favor of using entanglement as a diagnostic tool for 
integrability, independent support is needed to rule out 
possible bias due to the special relationship between mul- 
tipartite entanglement and the c-basis, or to the specific 
entanglement measure chosen. A natural option for cir- 
cumventing such limitations is to invoke GE. In partic- 
ular, the behavior of fermionic GE (as quantified by the 
-fu(L)-purity defined in Eq. P^ ) as a function of delocal- 
ization in the c-basis is depicted in Fig. [TO] for the same 
clean ID and 2D models examined above. Again, the 
observed delocalization dependence appears sensitively 
more pronounced in the 2D chaotic CcLSG, clS opposed to 
the ID integrable system. It is intriguing to think that 
this behavior ultimately reflects the fact that neighbor- 
ing energy eigenstates possess common features in the 
chaotic regime, as argued by Zelevinsky and coworkers 
[50| (see also t96i|). Although these findings reinforce 
the conjecture that, at least within the TBRE, a new 
entanglement-based signature of chaos may be identified 
in this way, validating this claim and its physical inter- 
pretation in more generality requires a dedicated investi- 
gation which we plan to report elsewhere. 



VII. CONCLUSION AND OUTLOOK 

We have provided a comprehensive quantitative anal- 
ysis of spectral properties, delocalization, and entangle- 
ment for the eigenvectors of disordered spin- 1/2 systems 
with Heisenberg interactions in one and two spatial di- 
mensions. Disorder in the system has been introduced 
through random applied magnetic fields, random inter- 
actions, or both. Our main findings and conclusions may 
be summarized as follows. 

(i) Although correspondence with non-integrability is 
well documented, the interpretation of LSI data is, as ex- 
pected, non-trivial in symmetry transition regions where 
invariant subspaces are partially overlapping - causing a 
tendency towards the Poisson distribution regardless of 
whether integrability is approached or not. 

(ii) Standard NPC measures for state delocalization de- 
pend entirely on which basis is chosen to represent the 
eigenvectors of the Hamiltonian, thus they need not de- 
tect the transition from integrability to chaos if the lat- 
ter is not accompanied by a significant change in the 
amount of delocalization. We have examined the effect 
of basis choice in two distinct bases of limiting integrable 
Hamiltonians as well as introduced a measure of relative, 
disorder-dependent delocalization which is shown to suc- 
cessfully detect quantum chaos in the presence of disor- 
der. We find that delocalization, particularly in the com- 
putational basis, does not achieve the level predicted by 
RMT even in chaotic regions for the Heisenberg model. 
Low connectivity and the localizing nature of the Ising 
pairing are likely to play a role in accounting for such 
disagreement. 

(iii) In delocalized systems, NPC and multipartite qubit 
entanglement show a dependence on energy that resem- 
bles the one observed for the density of states in the 
TBRE, refiecting a strong correlation among the three 
quantities. This complements findings on the infiuence of 
LDOS properties to dynamical aspects of entanglement 
in disordered qubit systems [23j . 

(iv) In the case where only random Heisenberg interac- 
tions are present, an interesting connection between level 
statistics, average eigenvector entanglement and delocal- 



14 



ization, and the presence of a random singlet phase was 
uncovered in ID. This may warrant further study, also in 
view of better elucidating what mutual relationships (if 
any) exist between quantum criticality and the ability of 
a system to sustain quantum chaos. 

(v) Delocalization and multipartite entanglement both 
show a substantial increase during the transition from a 
localized integrable Hamiltonian to a chaotic one. Simi- 
lar to NPC measures, both entanglement and GE measures 
are unavoidably relative to a choice of "local" resources 
- as captured by preferred subsystems or observables. 
However, comparison between two quantities sharing the 
same relative origin - NPC and P„ - results in a distinc- 
tive relationship which may serve as a quantum chaos 
indicator. Physically, such a relationship suggests that 
both entanglement and delocalization in a given local ba- 
sis essentially capture the same information about eigen- 
vector structure. Both the validity and implications of 
this potential entanglement signature, as well as a more 
substantial use of genuine GE, are points deserving addi- 
tional in-depth investigation. In particular, a promising 
venue to explore is the possibility that GE measures rel- 
ative to appropriate observable sets may allow to detect 
a transition to chaos starting from any integrable model, 
irrespective of its local or non-local nature with respect 
the the original operator language. 

Finally, from a technical standpoint, the general 
method suggested in and explicitly illustrated here 
for computing typical entanglement properties of random 
pure states localized to a proper subspace of states in 
Hilbert space is likely to find broader applications within 
both quantum chaos and QIS. 



APPENDIX A: TYPICAL SUBSPACE 
ENTANGLEMENT: EXPECTED VALUES OF 
BLOCK-PURITIES 

Given the general expression of relative purity in 
Eq. (|10p , each P„ may be expressed as a sum of squared 
expectation values over a normalized, orthonormal trace- 
less basis {Bi} for the distinguished observable set h. 
That is. 



where k is a normalization factor ensuring that 
max|^^{P„} = 1. Thus, the average purity over any en- 



semble of states is given by 



The method we shall follow to compute Pn will be based 
on the following result (proved in |93l|): 

Let Ti. be a finite- dimensional Hilbert space, with di- 
mension N , and let A be any traceless real symmetric op- 
erator on 7i, normalized such that tr{A'^) = N . Then for 



an ensemble of pure states \^) with real components taken 
uniformly with respect to the Haar measure on 0{N), the 
following relationship holds: 



N- 



For random states limited to a proper subspace, S C 
Ti., additional care must be taken because A need not 
be normalized or traceless after projection onto S. Let 
N'=dim{S) and let 11 denote the projector onto S. If we 
let n^n = aA'-\-l3t, where tr(^') = and tv{A'^) = N' , 
it follows that 



(nAn}2 



2a2 



N' 



0' 



The coefficients a and /3 may be determined from 
tr(nAn) and tr((nAn)2), 

2 tr(nAn)2 2 tr((nAn)2) tr(n^n)2 

p = , a = 



iV'2 



N' 



iV'2 



respectively. Therefore, in order to calculate P„, the 
trace norm of both the projection UBiH, and the pro- 
jection squared, (115^11)^, of each basis operator onto S 
must be determined. 

Since in our case each P„ is an average over the purities 
of {L/n) n-qubit subsystems, it suffices to determine the 
average subsystem purity of a single n-qubit subsystem. 
A convenient operator basis is provided by the set of n- 
qubit Pauli operators, that is, all products of single-qubit 
Pauli operators and the identity, Bi = aai ... (8) cri"', 
where ak — 0,x,y, or z. The choice of the identity, ctq, 
acting simultaneously on all qubits is excluded. It will be 
useful to consider the representation of the Pauli opera- 
tors as matrices expressed in the standard c-basis. The 
subspace S — Ho of states with total Sz = angular mo- 
mentum (no net magnetization) is relevant to our study, 
N' — Nq. The Pauli operators that have a non- vanishing 
projection into Ti.Q are: 

(1) those consisting of only cr^ and do operators; and 

(2) those consisting of an even number of ax and an 
even number of ay operators, with the remaining opera- 
tors being any combination of a^ and do- 

First, consider case (1). Each such operator is diagonal 
in the c-basis, hence it remains diagonal after projection 
into Ho- Because the eigenvalues of each such operator 
are ±1, ti{IlBiIl) may be determined by simply counting 
the number of c-basis states spanning Tip that correspond 
with each eigenvalue. For a Pauli string consisting of m 
az and (n^m) ao operators, the number of c-basis states 
with fc I's for the qubits acted on by the az operators is 

(T) 



im-k)- Thus, 



tr(ni?.n) = ^(-i)'=('^) 



L ^ m ^ 
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If m is odd, then the above quantity is 0. 

Next, consider case (2). Note that the effect of the 
operator ^ or ay ^ acting on a c-basis state is to flip 
the i'th qubit from to 1, or 1 to 0. Thus, for a c-basis 
state to remain in Ho after the action of a Pauh opera- 
tor, the combined number of ax and ay operators must be 
even. Additionally, since every state in the ensemble has 
only real components in the c-basis, those Pauli opera- 
tors that have an odd number of ay have zero expectation 
value for each state in the ensemble. For a Pauli string 
with a total of m ax and ay operators, m O's and I's 
will be flipped when the operator acts on a c-basis state. 
Thus, the qubits acted on must have an equal number 
of I's as O's in order to remain in Hq- There are (^2) 
possible assignments of O's and I's for such qubits. The 
remaining qubits (that is, those not acted upon by ax 
or ay operators) must also have an equal number of I's 
as O's since the state is in Ho- There are {^^j^S^/^} pos- 
sible assignments for such qubits. Therefore, there are 
(m/2) {(L-rn)/2) ^^'trix elements in such a Pauli opera- 
tor. Since each matrix element is 1, it follows that 



tr((nBin)2) 



m 
m/2. 



L — m 
(i-m)/2 



Since Bi acts non-trivially on every c-basis state, there 
arc no diagonal elements hence tr(ni?;n) = 0. 

In order to properly normalize purity on an n-qubit 
subsystem, consider the pure c-basis state |0, 0). Only 
Pauli operators consisting of a^ and a^ operators have 
non- vanishing expectation values for this state, and each 
such expectation value is 1. There are (2" — 1) such 
operators, therefore k = 1/(2" — 1) in this case. 

With the above ingredients, P„ may be calculated by 
determining how many Pauli basis states for an n-qubit 
subsystems fall into each category described above in re- 
gard to trace norm and squared trace norm upon projec- 
tion onto Ho . We consider various coarse-grained purities 
separately. 

• Pi: Only a'^^ has a non- vanishing expectation 

value. It has an odd number of az operators, thus 
tr((nc72n)2) = N and tr(ncr^n) = 0. Therefore, 



• P3: There are (^) Pauli operators containing m a^ 
and [n — m) uo operators, and in addition there are 
2'^"^^'' (2) PS'Uli operators containing 2 ax and {n — 2) 
(To and az operators and the same number of Pauli op- 
erators containing 2 ay and (n — 2) uo and cr^ operators. 
This yields: 



1 \Nq + 2 



7- 



3A| 



24 / L-2 
+ ]VoV(i-2)/2 



+ 



3A| 



• P4 : All of the classes of Pauli operators listed for P3 
must be considered for P4 also, with the following ad- 
ditional classes: the '^"^'^{22^1-4) P^^ili operators with 
exactly 2 ax and 2 ay operators; the 2"~^(2) Pauli oper- 
ators with exactly 4 ax operators; and the 2"~* (2) Pauli 
operators with exactly 4 ay operators. This yields: 



pr - 1 [ 2 
^ ~ 15\A/-o + 2 



15 



6A| 



48 



^ No 



L-4 
(L-4)/2 



6A| 
iV2 



Ni + No 

+ A 

N^ 



L-2 
{L - 2)12 



• Pq : The additional classes of Pauli operators to con- 
tribute to Pg are: the 2"~^(^ 2 n-e) P^'Uli operators with 
exactly ^ax and 2 ay operators; the 2"~^(2 4"_g) Pauli 
operators with exactly 2 ax and 4 ay operators; the 
2""^(g) Pauli operators with exactly 6 ax operators; and 
the 2"~^ (g) Pauli operators with exactly 6 ay operators. 



This yields: 



63 \. 



63 



15Ai 



+ 



A^o + 2 

480 / L-2 \ 480 



15A| 



Nl 



No 



32 / 

nMl 



{L-2)/2) 
L-6 
6)/2 



No 



L-4 
(i-4)/2 
15A| 
Nl 



3. 

Nl 



Letting L = 12, No = 924 in the above formulas leads to 
the n-local purity values quoted in Table I. 



Pi = 



No- 
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